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Abstract 

Long-time asymptotic of field-field correlator for radiation propagated through 
a medium composed of random point-like scatterers is studied using Bete-Salpeter 
equation. It is shown that for plane source the fluctuation intensity (zero spatial 
moment of the correlator) obeys a power-logarithmic stretched exponential decay 
law, the exponent and preexponent being dependent on the scattering angle. Spatial 
center of gravity and dispersion of the correlator (normalized first and second spatial 
moments, respectively) prove to weakly diverge as time tends to infinity. A spin 
analogy of this problem is discussed. 
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1. Introduction 



The fluctuation phenomena in scattering of electromagnetic waves propagating in ran- 
dom media was the subject of extensive studies during the last decade. This activity 
was inspired by experimental evidence for a number of the phenomena (e.g. enhanced 
back-scattering [1,2], spatial and frequency correlation of back-scattered and transmit- 
ted radiation [3,4]; for a review see [5]) as well as by realizing a close similarity of these 
phenomena to fluctuation phenomena in electron transport in disordered systems (meso- 
scopics) [6,7]. The latter brought new physical insight and theoretical methods into the 
field established in optics many years ago [8]. The study of fluctuation effects was ex- 
tended also to dynamic disorder in which the scatterers move at random. Golubentsev 
showed that the contribution of interference back-scattering to static field auto-correlator 
(in particular, average intensity), which is essential for immobile scatterers, is suppressed 
when they move [9]. The similar result for temporal-spatial auto-correlator (dynamic 
coherence function) was obtained by Stephen [10]. Short-time fluctuations were studied 
using diffusion approximation in [10] for the cases of a point source in an infinite medium 
and a plane wave incident on semi- infinite medium (Milne problem). However, the valid- 
ity of this approximation for time-dependent disorder remained questionable and deserves 
special study, to which the present paper is devoted. 

2. Definitions and Statement of the Problem 

Let the electric field E (r, t) of light obeys scalar wave equation with a source cur- 
rent density j(r,t). Coherent field (CF) and field coherence function (FCF) are defined 
respectively by 

E c (r,t) = {E(r,t)), (1) 

T (n, h; r 2 , t 2 ) = (E (n, h) E* (r 2 , f 2 )) - (E (n, h)) (E* (r 2 , t 2 )), (2) 

where brackets mean average over realizations of underlying random process. In the 
case under consideration it is the scatterers motion that causes the randomness. Point- 
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like scatterers are considered, so individual scattering can be described by scattering 
amplitude f s (see e.g. [8], vol.1). In weak multiple scattering limit A <C I (where A is the 
radiation wavelength, I = 4:7i/n s \f s \ 2 is radiation mean free path, and n s is the scatterers 
density) CF is easily calculated. For space-time Fourier transform of CF one obtains 

E c (p, w) = (p, u) G (p, u) , G(p,u) = — -I -, (3) 

c p 2 - k£ - n s f s 

where k = 2n/\ is the wave number, oj = ck is the radiation frequency and G(p,cu) 

is 'one-particle' Green function (see e.g. [6,7]). When considering FCF slow coordinate 

r = (ri + r 2 )/2 and time t — (t± + t 2 )/2 as well as fast ones: p — (i"i — r 2 )/2 and 

t = (ti —t 2 )/2 are introduced. Then Fourier transform of FCF (Eq.2) with respect to the 

fast variables 

- . . f f°° . . ^ . . d 3 pdr . A . 

r(r,*;p,w) = / / exp(ttJT-zp-p)r(ri,fi;r 2 ,t 2 )— — 4 (4) 
J J-oo 

is commonly used. This quantity may be interpreted as energy density in incoherent 
portion of the radiation with the frequency u and propagation direction n = p/p at the 
point r and time instant t (see e.g. [8], Vol. 2). 

Main issue of the present paper is the study of f (r, t; p, a;) at long times neglecting the 
coherence effects. By 'long times' I mean the regime when A (t), mean square displacement 
a scatterer moves during time t, is lagrer than the wavelength: 

A (t) > A. (5) 

3. Bete-Salpeter Equation 

Bete-Salpeter equation in weak individual scattering regime was derived in [10]. In 
general case, except some parametric factors, the structure of this equation remains un- 
changed. For f (q, t; p), spatial Fourier transform of f (r, t; p), it reads 

f (q, t; p) = n s \f s \ 2 G (p - q/2) G* (p - q/2) x 

d 3 p 1 



exp 



-|A(t) 2 |p- Pl | 2 



r (q, t; Pl ) - E c ( Pl - q/2) E* c ( Pl + q/2) . (6) 

(2tt) 



The argument u> here and further is omitted for brevity. It is worth noting that Eq.6 may 
be also derived using Foldy-Twersky ansatz (see [8], Vol. 2). At distances much larger 
than wavelength, that is at q/k <C 1, the approximation 

irl 8(\p\-k ) 



G (p- q /2) G .(p- q /2)- 2 , o2 ..^.^ 
holds well, and Eq.6 with account of Eq.3 will be satisfied by the substitution 

f (q, t; p) = (|p| - fc ) 7 (q, t; p/k ) 



(7) 



(8) 



that transforms Eq.6 to the transport equation for the on-shell FCF 7 (q, t; n) 



(1 - iZq • n) 7 (q, t; n) - / exp 

J\m\=l 



^1 I 

n — rii 
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7(q,*; n i 



47T 



exp 



|n |=l 



^1 I 

n — n 
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j (/c n - q/2) j* (k n + q/2) dfl 



1 — i/q • n 



47T 



, J=(k A) 2 . (9) 



For point source Eq.9 coincides with that derived by Golubentsev [9] using a special 
diagram technique whose rules take explicitly into account Eq.7. The sources of other 
types were considered [11] within real-space approach (RSA) after Rosenbluh et al. [12]. In 
the next section a connection between Bete-Salpeter equation and RSA will be discussed. 

4. Real-Space Trajectory Decomposition 
and Spin Analogy 

Solving Eq.9 by iterations of the integral term represents its solution as an average of 
7 (q, t; n , n), source-detector directional FCF, over the propagation directions n allowed 
for radiation 

j (k n - q/2) j* (k n + q/2) 7 (q, t; n , n) (10) 

|n |=l 

with the series expansion in the integrand: 

7(q,t; n ,n) = 
4 




(47T) 



,N-l 



N-l ' 



11^ = 11. (11) 



Here Wn are scattering configuration weights given by 



N 



1 



1 1 A (q; Hn. 11, H \ ) = J J 

3=0 



1 — iZq • n, ' 



(12) 



At iV 3> 1 and |r| 3> / the inverse transform Fourier of this product is strongly peaked at 



trajectories terminated at the point r. Each jth (0 < j < N) straight-line part of them 
has the length I and directed along the vector n,. The initial n and final directions 
are determined by a source and a detector, respectively, and each change of 'free-flight' 
direction occurs after a scattering act. Scattered directions ni, . . . , n^-i are distributed 
with correlation if J ^ 0. 

Thus, basic assumptions of RSA [12] result from the Bete-Salpeter equation. In RSA, 
however, Wn are calculated using apriori random-walk considerations, and the exponen- 
tial weights due to scatterers motion are approximately averaged by averaging of their 
exponents [12]. An analysis shows that W N of RSA is an average of right-hand side of 
Eq.12 at ql <C 1 and N ^ 1 over fully random n^. With the only this less stringent ap- 
proximation 7 (q, t; n , n) will coincide with conditional partition function in an ensemble 
of one-dimensional classical Heisenberg spin systems, each with the same exchange inte- 
gral J ', but different numbers iV > 2 of species and weights Wn- This analogy facilitated 
an estimation of FCF within RSA making use the methods of magnetism [13]. 

5. Long-Time Asymptotic of Coherence Function 

Scatterers motion removes diffusion singularity: 7(q, t; n , n) oc q~ 2 at q — > specific 
for immobile scatterers. At ql 1 it makes sense to consider hydrodynamic expansion 



r = Z£ 




rij, so in real space FCF appears to be a sum of path integrals over polygon 



■oc 



7 (q, t; n , n) = # (n , n) + i k 




fi'ai...a fc ( n 0; n ), 



(13) 



fc=l a 
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where a's are the Cartesian indices, and the functions g ai ...a k (*! n o, n ) 

^<9 fe 7 (q,*l n ,n) 
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a\...a k 



(n ,n) = % 



dq ai ...dq ai 



q=0 



= y r Q1 ...r ttfc 7(r,t; n ,n) d 3 



(14) 



are spatial moments. The first term g ( n o, n ) of the expansion in Eq.13 (the zero moment) 
is the intensity of the fluctuations in a macroscopic volume. Next two moments normalized 
to the zero one have clear physical meanings. The vector and matrix 



9a (no, nj 



r a rfs 



g a p (n , nj 



(15) 



g {n ,n) g (no, nj 

are the spatial center of gravity and dispersion matrix, respectively. Eq.13 can be directly 
obtained by iterations of the term oc qn in the transport equation for 7(q, t; n , n), which 
differ from Eq.9 only in the right-hand side. Such a procedure creates recurrence equations 
for the moments, which enable us to calculate them recursively, but analyzing the full set 
of the moments (same as solving the original equation) is hardly possible. Nevertheless, 
the asymptotics of the above mentioned three moments at J — > oo can be calculated 
rigorously. This limit is quite reasonable since Eq.5 implies the strong inequality J >> 1. 

5.1. Fluctuations Intensity 

Using Eq.10 and Eq.9 one obtains the equation for the zero moment 
(1 - £)go (n ,n) = exp 



J \ 12 

n — n 

2 1 01 



(16) 



where £ is the same integral operator as appears in the left-hand side of Eq.9. The eigen- 
functions of the operator £ is well known spherical harmonic (n) with the eigenvalue 
Xl degenerated with respect to the azimuthal number m = 0, ±1, L: 



£Fi m) (n)= J exp 



n - nJ 

2 1 11 



' ( ""'«* = A# (n). 

4"7T 



rr (ni) 



(17) 



|m|=l 



Applying the addition theorem for spherical harmonics and the eigenfunction expansion 
of the solution to Eq.16 one gets closed expressions for the eigenvalues: 



X L = A[L + -,J), A(v,x) 



71 

2x 



exp (-x) I u (x) , 



(18) 



and the zero moment: 

g (n , n) = J2( 2L + ^Y^T Pl {C0S9) ' (19) 

L=0 L 

where I v (x) is modified Bessel function [14], Pl(x) are the Legendre polynomials and 9 is 
the scattering angle: cos 9 = n ■ n . 

Apparent fluctuation intensity depend on the detection direction n if the source is 
directional (i.e. a narrow range of the radiation directions n is allowable). In the opposite 
case of point source all terms except isotropic one with L = are averaged out when 
substituting Eq.19 into Eq.10 at q — 0. It is the value of J that determines how strong 
the anisotropy of Eq.19 is. Power expansion of I v {x) with respect to x leads to Al = 
J L /(2L + 1)!! + O (J L+2 ) at J < 1. Thus, at short times Eq.19 is dominated by 
isotropic term oc J7~ 2 , and the anisotropic terms are relatively small as J L+2 . By an 
analogy with the case J = 0, q ^ this regime may be called diffusional in the sense 
that multiple scattering retains the isotropy of individual scattering on macroscopic scale. 
Quite opposite behavior takes place at long times. Using the asymptotic of I v (x) at x ^> 1 
one obtains \l = (2J") 1 + O {J~ 2 ) at J 3> 1. Independence of the leading term in Al 
on L makes all the terms in the series (19) important. Moreover, since next-to-leading 
asymptote of Al contains the factor L (L + 1) the terms with large orbital numbers seem to 
be even more important. This means that the calculation of go(n , n) requires summation 
of the whole infinite series. As a result, scattering angle dependence will be strong. 

To exploit classical summation technique the series of interest is represented as a 
residue sum over points Z L = i (L + |), the zeros of cosh^ 

oo 

g (n ,n) = 2ni ^ Res 

L=0 

where •& = n — 9, and P u (x) is the spherical function. Adding the sum of residues at z h 
the zeros of the function 1 — A (—iz, J) in upper half-plane, to Eq.20 gives the integral 



— izA (—iz, J) 
l — K(—iz,J) cosh7rz 



P_i 

2 



(cosi?) 



(20) 



z = 7,t 
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of the function in square brackets over the line (—00 + i ■ 0, +00 + i ■ 0), and therefore 
, f +QO -i(A(-i(,J) P-i-g(cosfl) *iP-*-*(cos0) 



l-A(-i(,J) coshvrC "' ^ cosh{n Zl )f z A {-iz t , J) 

Here A (— i(, J) for real C means the value on the upper branch of real-axis cut. Because 



A (—iz) — A (i z) both zi and — zi are the roots. Using this and the properties of the cone 
function V _i_ iz (x) [15] one obtains after a transformation 



+°° CImA (-<, J) P_i_, c (cos^) 

[1 -ReA(-t(,J)} 2 + [ImAHCjT)] 2 ' coshvrC 



EziP_i_ iz . (cos$) 
. (21) 
^ >0 cosh (t^A (-**,, J) 

With this exact representation the desired asymptotic can be now handled. Eq.18 and 
Bessel functions theory [14] give at Imz > 

7T 

A (— iz, J) = — /=== J ex P [— J" (1 — cos-u)] cosh (uz) du 



00 

.sinh(vr,2) f r „ . , , 

+ i — / exp [-J (1 + cosh w) + iwz] dw. (22) 

\/2-kJ J 


Both integrals in Eq.22 for real C can be estimated at J — > 00 using the Laplace method. 
At [CI > v/2^7 one obtains 



6XP ' 2J J + ^ sinh «) ex P ( ~ 2 ^~ 2 7 



(23) 



Let us assume, in addition, that |£| <C 2J/n. The integral in Eq.21 is now easily calculated 
because under this condition lorentzian-like factor in the integrand has exponentially small 
broadening ImA (—i(, J\ and hence is well approximated by the delta-function 

^ (C - Co) 

IfReAHCcJ)!' 

where Co is the positive root of real equation ReA (— i(, J) = 1. Substituting here the 
real part of Eq.23 gives the asymptotic of the root 

Co ^ y/2J\og{2J). (24) 



This equation is consistent with the assumed inequalities provided that 

V2J 



7T 



> y/log(2J) > 1. (25) 



An analysis based Eq.22 shows that in this limit the equation A (—iz, J) = 1 has in the 

upper right quadrant only one root z given asymptotically by 

.exp(7rCo-2J r ) 
z o-(o + * • (26) 

Due to the inequalities (25) the imaginary part of this root is exponentially small, and can 
be discarded when calculating the residue term in Eq.21. Then the integral and residue 
contributions give together 

P_i_, Co (costf) 

9o (n , n) ~ 6ttJ 2 <oV \ (27) 
cosh 7l(o 

At the end, using in Eq.27 the asymptotic of cone function at large argument [15], one 
obtains for hemisphere < 9 < it 

3o(n , n) ~ Unjf (9) exp [-9^/2 J log (2J)] , (28) 

where / (9) = 1 at backward scattering (9 = n), otherwise at 2ir sin 9 3> [2 J log (2J)]~ l l 2 

f{9) ~ (27rsin^)- 1/2 [2 J \og(2J)r l,A . (29) 

At forward scattering (9 = 0) the whole derivation is inapplicable because all terms of 
the series in Eq.19 are positive. In this case the fluctuation intensity proves to be greater 
than unity at any time. 

To express Eqs.28, 29 via t one may use the models of the scatterer motion described in 
the literature. In two simplest, but practically important, cases of ballistic and diffusive 
motion: 2J7" = (t/r a ) a , where a = 2 and a — 1, respectively. The formulas for time 
constants r a are found, for example, in [9,10]. Thus, at t — > oo the fluctuations decay 
according to power-logarithmic stretched exponential law 



#o(n ,n) ~ exp 



(30) 



5.2. The First and Second Spatial Moments 

Calculating the asymptotics for g a (n ,n) and g a p (n ,n) is very involved task. To 
make things easier one might introduce local coordinate system with z - axis parallel 
to vector n and consider local longitudinal components g z (n , n) and g zz (n , n) . Other 
components are zero if spatial variation occurs only in the radiation direction. This is the 
case for an infinite source concentrated at the plane z — 0. CF calculated using Eq.3 is 
one half of the sum of the two attenuated forward and backward plane waves 

/ , \z\ 



^(r)ocexp^-^, (31) 

that is there are two allowed radiation directions n = iIq = (0, 0,±1). Fortunately, the 
superposition principle holds also for FCF (Eq.10). This turns out to be equal to the 
sum of half-weighted 7 (r, £;n^,n), both being independent on x,y. Therefore, consid- 
ering CF and FCF for one direction, say n^, one could hope to emulate a transmission 
problem where plane wave impinges on the boundary z — 0. Both CF and FCF vanish 
in the counterpart half-space, and coherent portion of radiation propagates only forward, 
whereas the fluctuations develop in all directions with the intensity at t — > 00 given by 
Eqs.28, 29. 

The first longitudinal moment, calculated recursively as outlined above, is given by 



9z ( n o, n) = I g (n , n) + I 



L=0 



(L+1)A L+1 + L\l-i 



P ^"» 9 \ (32) 

1 — Ar 



1 — Al + i 1 — Al_i 

At 9 > and J — > 00 the summation technique developed above can be applied to Eq.32 
as well. After division by Eq.27 one obtains the asymptotic for normalized first moment 



z~^ j = J 2J - (33) 
" 6 ' I \j\o g (2J) [alog(t/r a )] 1/a 

tending to infinity if < 9 < ir, and ~z~ — > if 9 = n. Closed formula for the second 

longitudinal moment g zz (n , n) is derived along the same lines as Eq.32, and long-time 

asymptotic of normalized second moment is obtained making use the summation method 
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developed above. Expression for g zz (n , n) is rather lengthy, and obtaining the asymptotic 
involves processing many cumbersome terms without bringing new ideas. This analysis 
is omitted, and the final result is as follows 

T ^ejl. mi ( 34) 

" 3 3 alog((/T„)' [<> > 

In conclusion, whereas the fluctuations intensity decays at t — > oo according to Eq.30, 
spatial center of the FCF 'runs away' according to Eq.33, and due to Eqs.33, 34 spatial 
dispersion of the fluctuations 5 = \J z 2 — (z f ~ £ (1 -sm 2 9/12) 1/2 /V3 also tends to 
infinity but with much weaker anisotropy. In other words, at long times random motion 
of the scatterers weakens the fluctuations, thus tending the radiation to a coherent state, 
but enhances drastically the spatial range of the field correlations £ over its value for static 
disorder ~ I. 

Acknowledgements. 

I am grateful to Dr. E. Kogan, who attracted my attention to this problem, for 
enlightening discussions and critical reading of the manuscript. 



11 



References 



[I] M.P. van Albada and A. Lagendijk, Phys.Rev.Lett., 55, 2692 (1985) 
[2] E. Wolf and P. Maret, Phys.Rev.Lett., 55, 2696 (1985) 

[3] A.Z. Genack, N. Garcia and W. Polkosnik, Phys.Rev.Lett., 65, 212 (1990) 

[4] M.P. van Albada, J.F. de Boer and A. Lagendijk, Phys.Rev.Lett., 64, 2787 (1990) 

[5] M. Kaveh, Waves in Random Media, 1, S121 (1991) 

[6] B. Shapiro, Phys.Rev.Lett. 57, 2168 (1986) 

[7] M. Stephen and G. Cwilich, Phys. Rev. B 34, 7564 (1986) 

[8] A. Ishimaru, Wave Propagation and Scattering in Random Media, Vol. 1-2 

(Academic, New York 1978) 

[9] A. A. Golubentsev, Sov. Phys. JETP, 59, 26 (1984) 
[10] M. Stephen, Phys. Rev. B 37, 1 (1988) 

[II] I. Edrei and M. Kaveh, J. Phys. C 21, L971 (1988) 

[12] M. Rosenblue, N. Hoshen, I. Freund and M. Kaveh, Phys.Rev.Lett., 58, 2754 (1987) 

[13] E. Kogan, private communication; 

M. Auslender, Bull. Israel Phys.Soc. 39, 134 (1993) 

[14] G.N. Watson, A Treatise on the Theory of Bessel Functions, 

(University Press, Cambridge 1966) 

[15] V.V. Lebedev, Special Functions and Their Applications, 

Selected Russian Publications in the Mathematical Sciences, ed. by R.A. Silverman 
(Prentice-Hall Inc, Englewood Cliffs, N.J. 1965). 



12 



